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SUMMARY 


Radiation- interchange conf iguration factors are derived for axisym- 
metrical sections of cylinders, cones, and hemispheres radiating inter- 
nally to annular and circular sections of their bases and to other axisym- 
metrical sections. The general procedure of obtaining configuration 
factors is outlined and the results are presented in the form of equations, 
tables, and figures. 


INTRODUCTION 


The high temperatures which are encountered in components of aero- 
space vehicles have brought about a renewed interest in heat transfer 
by radiation. For certain components, the heat transferred by radiation 
can be shown to overshadow that transferred by conduction. In the present 
paper axisymmetric radiation for some common axisymmetric shapes is 
studied. 

Radiative transfer of heat from one area to another depends, among 
other things, upon the fraction of the radiant energy emitted by one area 
which is intercepted by a second area. This fraction is identified by 
several names, such as the configuration factor, the interchange factor, 
the angle factor, or the geometric view factor, and is a function of the 
geometrical relation of the areas involved. In the present paper, the 
term configuration factor will be used to designate the fraction. 

Configuration factors are available for radiation between various 
surfaces (see refs. 1, 2, and 3 ) "but, for the most part, the areas which 
are involved are plane. The purpose of the present paper is to provide 
configuration factors for some of the more common nonplanar surfaces. 

Some of the configuration factors presented herein are obtained, in 
appendix A, by integrating the basic equation which defines the factor 
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and the remainder are obtained, in appendix 3, by the application of 
configuration- factor algebra. The technique!', of configuration- factor 
algebra make it possible, in some situations^ to obtain the desired 
configuration factor from available factors without the need for 
integration. 

In addition to being listed in tables and given in the form of 
equations, the results are presented in carpet plots which permit an 
estimate of the magnitude of a given factor end show the effect of 
varying the proportions of the surfaces involved. 


SYMBOLS 


A area 

a radius of the base of a surface of revolution 

C circular area 

F configuration factor defined by equation (4) 

H height of a cone 

j,k,m, n integers 

4 length between the mth and nth places 

M nondimens ional parameter, r^ya 

N nondimens ional parameter, L^/a 

q energy per unit time 

area of ring between the mth and nth planes 
r radius of circle 

5 distance between centers of the areas exchanging radiant 

energy 

T absolute temperature 

x,y, z Cartesian coordinates 


hI ONO\OJ 
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p,0,z polar coordinates 

p,0,(p spherical coordinates 

a Stefan- Boltzmann constant 

half the apex angle of a cone 

i|r n angle between the normal to the area An and the line between 

centers of the area An and the area which intercepts 
radiation from A n 

0 ) solid angle 

Subscripts: 


0, 1, 2, 3, j,k, n identification of an area, plane, or point 
Cj circular area in the base of a body of revolution 

C 




an annular area in the base of a body of revolution 

( C J>k = C j " C k) 


dApcLAg from an area cLA-^ to an area CLA 2 

j,k from an area J to an area k 

Superscripts: 

1,2, 3, a,H,m, n identification of an area, plane, or point 


ANALYSIS 


Black- Body Radiation Between Two Isothermal 

Surfaces of Arbitrary Orientation 

Consider the exchange of radiant energy between two isothermal 
black surfaces, A^ and A 2 , of arbitrary orientation separated by a 
nonabsorbing medium as shown in figure 1. The energy per unit time 
leaving the first elementary surface dA^ in the direction of the 

second dA 2 is given by (see ref. l) 

d< JdA 1 = i *1 cos +1 dAi 


( 1 ) 
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where 


q energy, per unit time 

O Stefan- Boltzmann constant 

TjL absolute temperature of dA 1 

angle between normal to dA^ and line between centers of 
areas dA-^ and 

The portion of the energy per unit time leaving dA^ which is 
intercepted by dA 2 depends upon the solid ingle dco subtended by dA 2 
and can be expressed as 

dq dA 1 ,dA 2 = * T 1 cos h (2) 


where 


do> = 


cos i|»2 


cLAg 


and S is the distance between the centers cf the areas dA-j_ and dAg. 

The energy per unit time which leaves the surface dA-^ and is 
intercepted by the surface dA2 (eq. (2)) can therefore be expressed as 


. cos ^ cos + 2 
' ? T 1 p “l ^ 


( 5 ) 


By defining 


A 1 ’^ " * A i Ja 1 Jaz 


cos i|r^ cos 


dAg dA-^ 




equation (3) becomes 


%A2 = aT ^ F 


ITApAg 


( 5 ) 


no vo vo tr* 
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The term 


F A p A2 


depends only upon the geometric configuration of the 


two surfaces and is known as the configuration factor. The configura- 
tion factor is defined as the fraction of the total energy per unit time 
which leaves a surface and is intercepted by a second surface. The above 
procedure can be repeated to determine the energy per unit time which 
leaves the second surface and is intercepted by the first with the fol- 
lowing results: 


%,A 1 ' aT 2 A 2 F A 2 ,A 1 


( 6 ) 


where 


■A^A-, 


1 

it 


cos 4^ cos ^ 


dAi (IA2 


It can be seen from equations (4) and ( 7 ) that 


At F 


= A^F„ 


1 t A 1 ,A 2 "2 i A 2 ,A ] 


(7) 


(8) 


Equation (8) is known as the reciprocity theorem and, as is seen 
in appendix B, is very useful in the application of configuration- factor 
algebra. For brevity, whenever the areas involved are understood, equa- 
tion (8) is written as 


A 1 F 1, 2 “ A 2 F 2, 1 

The net exchange of radiant energy between Aj^ and of fig- 

ure 1, obtained as the difference between equations (5) and (6) and 
simplified by the reciprocity theorem, is 

Vt * 0A l F l,2( T i - T 2) (9) 


Black-Body Radiation in Closed Systems 

Consider a closed system composed of n isothermal black-body 
surfaces separated by a nonabsorbing medium. The net heat flow result 
as presented for two isothermal black-body surfaces can be extended to 
include the n isothermal black surfaces in the following manner. 

According to the Stefan-Boltzmann law, the radiant energy leaving 
an isothermal surface j is 
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The radiant energy incident upon the area from all other surfaces 

in the system is 


I j 

k=l 


(ID 


so that the net exchange of energy becomes 


Inet 


- « t j a j 


n 

I 


K A 


k=l 


( 12 ) 


Finally, equation (12) can be reduced by applying the reciprocity theorem 
to obtain the following equation: 


I 


^net 


a A A 


n 


-i 

lr=l 


k J,k 


(13) 


With a knowledge of the configuration factors F. , , equation (13) 

J> & 

can be used to obtain heat flows or temperature distributions in a closed 
system. Reference b demonstrates the use of equation (13) when radiant 
heat transfer is accompanied by heat transfer by conduction. 


Configuration Factors 

In practice, configuration factors can be obtained by experimental, 
numerical, and analytical means. (See, for example, ref. 1.) In the 
present paper, some configuration factors are found directly from equa- 
tion (4) and some indirectly from equation ( 4 ) through configuration- 
factor algebra. 

The shapes considered are surfaces of revolution (cylinder, cone, 
and hemisphere) with the ends closed by plane surfaces. All areas con- 
sidered are axisymmetrical and therefore the resulting configuration 
factors are applicable only to surfaces exhibiting axisymmetrical tem- 
perature distributions. 

In the study of radiant heat transmission within a system composed 
of a body of revolution and a base plane six types of general configuration 
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factors are required. These six types of configuration factors, which 
are derived in the present paper, are shown in figure 2 for the case 
where the surface of revolution is a cylinder. Similar configuration 
factors have been determined for the cone and hemisphere. 

In all six configurations, the surface of revolution is divided into 
rings by one or more planes which are parallel to the base at heights of 

1^, L^, . . . l£. The subscripts and superscripts indicate the planes 

between which the length is measured (zero being the base plane) so that 

12 n 

the rings between these planes will be known as Rq> Rq, • • • Rq- 

Areas in the base plane considered are either circular or annular 
and are designated by C q, C-^, ... C n and by S n -±,n> res P ec ‘ki ve ly* 

The circular area Cq represents the full base of the cylinder so that 
C ^ > C n . The annular region C n-1 ^ n represents the region between the 
circular areas C n _-^ &nd C n so that 

c n-l,n - Vl - C n 


Since the circular and annular areas are normally in the base plane, 
there is usually no need to specify the plane in which they lie. How- 
ever, for a few cases it is necessary to specify the plane and this will 

be done by the use of superscripts. For example, = ^n-1 “ ^n 

so that the annular region is in the mth plane and is equal to the area 
contained between concentric circles in that plane. Whenever C terms 
appear without superscripts the area is understood to be in the base plane. 

By using this method, Rq, would indicate a ring on the surface 

of revolution extending from the base plane to the first plane above the 
base exchanging radiant energy with a circular area in the plane of 

the base. In the same manner R^, C- L g would indicate a ring lying 

between the first and second planes above the base exchanging radiant 
energy with an annular area, Cj_ - in the base plane. 


By using this symbolism, the six configurations presented in fig- 
ure 2 are designated as R^, C-^, RpCp 2' ®i^i 2^ ^2^0* 

p P 
and R^, R^. 
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The derivation of the configuration factors is presented in the 
appendixes. The configuration factors for r,he geometry of R 1 , c are 
obtained by integration of equation (4) for cylinders, cones, and hemi- 
spheres in appendix A. With the exception of r|,R^ for the hemisphere, 
the remaining configuration factors are obtained by using configuration- 


factor algebra and the equations derived for the geometry of R^, c . 

The use of configuration-factor algebra is explained and demonstrated 
in appendix B where configuration factors are given for the geometries 


of 




r 1 r 

V L 1,2' 


R 1 ,C 1,2' 


R 2> R q> aad 


2 2 

r i' r i- 


The geometry 


2 0 


for the hemisphere is not amenable to the use of configuration-factor 
algebra. This situation results from the fact that configuration- factor 
algebra depends to a large extent on dealing with similar surfaces and 

the geometry of R^jRq results in spherical segments which are not 

hemispheres. Therefore, the result for this case is obtained by integra- 
tion in appendix A. 


RESULTS AND DISCUSSION 


Since all the configuration factors given in the present paper are 

obtained with one exception (r|Rq for a hemisphere j from three primary 

equations (A7), (Al6), and (A26) involving tae geometry of R^^ and 

configuration- factor algebra, only evaluations of these three equations 
will be discussed in any detail in this section. The configuration 

3 1 

factor derived in appendix A for R^Rq when the surface of revolution 
is a hemisphere is elementary and will not require discussion. 

Table I is a summary which indicates, by number, the proper equa- 
tion to use for the cases previously described. In addition to the 
specific surfaces of revolution treated in detail, fundamental equations 
are indicated for use with an arbitrary surface of revolution. 


Cylinders 

The configuration factor for the geometiy RqjC-^ when the surface 

of revolution is a cylinder is derived in appendix A and is given by 
equation (A7). Equation (A7) is given in a r.ondimensional form by 
equation (A8) which is 


^ o\ ON OJ 
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^ + 2N 2 (l + M 2 ) + (l - M 2 ) 2 - (l - M 2 ) - N 2 

1 

The nondimens ional parameters are 

M = r^j a 

B - L^/a 

where a is the radius of the base, r^_ is the radius of C^, and 
is the height of Rq. (See fig. 3-) 



The term M is a ratio of the radius of area C-j_ to the radius 
of the base of the cylinder whereas the term N is a slenderness 

ratio. 

Table II presents results from the nondimens ional equation (A8) 
for a wide range of r^J a and L^/a. The data of table II are also 
given in the form of a carpet plot in figure 4. 


Cones 

The configuration factor for the geometry of Rq, when the sur- 
face of revolution is a cone is derived in appendix A and is given by 
equation (Al6), which is a lengthy equation that results from the evalua- 
tion of a nonelementary integral. Because a large number of terms in 
the equation must be defined, it will not be repeated in the text and 
reference should be made to appendix A. Table III gives the results of 
the evaluation of a nondimensional form of equation (Al6) for combina- 
tions of I^H and r^|a between 0.1 and 1.0 and for cone half-angles 

of 5°, 10°, and 20°. The dimensions Lq and r^ are as shown in 
figure 5- 

The data of table III are presented in the form of a carpet plot 
in figure 6. Figure 6 contains three parts, one for each of the half- 
angles considered. As expected, the evaluation of equation (A17)> shows 
that the cone results approach the results obtained for the cylinder as 
the base angle approaches rt/2. 
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Hemispheres 

The configuration factor for the geometry of R^,^ when the sur- 
face of revolution is a hemisphere is derive! in appendix A and is given 
by equation (A26). Equation (A 26) is given in nondimens ional form by 
equation (A 27 ), which is 



1 



1 - M 2 ) 2 + l+M 2 !! 2 - (l - 



where the nondimens ional parameters are, as for the cylinder, 

M = r-^a 

H - L^/a 

_1 

"q, r.^, and a are as sshown in figure 7 . Table IV 

the evaluation of equation (A 27 ) for combinations 
between 0.1 and 1.0. 


and the dimensions 
presents results of 
of r-j^a and I^/a 


The data of table IV are presented in the form of a carpet plot in 
figure 8. Figure 8 shows that F_i „ is constant for all values of 

c : 

N when M = 1. 


CONCLUDING REMARKS 


Configuration factors are presented which can be used in heat- 
transfer studies involving nonplanar surfaces at high temperatures. 
These configuration factors have been derived for axisymmetrical sec- 
tions of cylinders, cones, and hemispheres rediating internally to 
circular and annular regions of their bases cr to other avl sy mme tric 
sections. Some of the factors were obtained by integrating fun dame ntal 
equations expressed in terms of convenient coordinates. The remainder 
of the factors were obtained by utilizing configuration-factor algebra 
and the results of the integrations. The use of configuration— factor 
algebra is explained and demonstrated. The calculated radiation con- 
figuration factors are given in tables and plots. 


Langley Research Center, 

National Aeronautics and Space Administration, 
Langley Air Force Base, Va. , July 20 , 1961. 
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APPENDIX A 

DERIVATION OF CONFIGURATION FACTORS FOR CYLINDERS, 
CONES, AND HEMISPHERES BY INTEGRATION 


The equation for the portion of the total radiation from an elemental 
area dAq which is incident upon an elemental area dA 2 is derived in 

the body of the paper as 


Y A 1 ,A 2 jtA-l J f 


cos cos +2 


Ag 


dAg 


dA-, 


(Al) 


where 

S distance between dA-^ and &A 2 


\|r^, ijfg angles between the line joining the areas dA]_ and dAg and 
the normals to the respective areas. 


The configuration factor will be derived for the most general case 
of a section of a body of revolution extending from the base plane to a 
given plane above it exchanging radiant energy with an axisymmetrical, 
circular region located in the plane of the base as shown in figures 3, 

5, and 7 for cylinders, cones, and hemispheres, respectively. The param- 
eters of equation (Al), ^ 011(1 s > must ^ expressed in terms of 

the chosen coordinates so that the equation can be integrated. 

In order to keep the solution as general as possible, the integra- 
tion of equation (Al) will be carried out over surfaces designated as 
A]_ and A 2 . The area A^ will represent the area on the surface of 

revolution R^ and the area A 2 will represent the circular area C- L 
in the base. This will apply throughout the derivation of the configura- 
tion factor for the geometry of Rq, for cylinders, cones, and hemi- 

3 1 

spheres. For the geometry of R^Rq for a hemsiphere, the areas A 1 
and Ag represent the upper and lower rings, respectively. However, in 
al 1 cases the limits will be written in general terms. 
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Configurations Factors for the Geometry Rq, C-^ 

Cylinder . - When the surface of revolution is a cylinder (see fig. 3), 
the parameters of equation (Al) in terms of th; polar coordinates p, 9, 
and z become 

a - p p cos/'e, - 9 ? \ 

cos t-L = ^ L (A 2) 

s 

cos if 2 = | (A3) 

S 2 = z 2 + a 2 + p 2 2 - 2ap 2 cos(9 1 - 9 2 ) (A4) 


dA-^ = a d9-^ dz 


dAp = P 2 dp 2 d0 2 


Integration over 9 2 , from 0 to 2n (after substituting eqs. (A2), 
(A3), and (A4) into eq. ?Al)), gives the configuration factor from the 
area dAq to the differential ring in the base 2np 2 dp 2 as 


dAiF 


1 dA 1 ,2«p 2 dp 2 


_ 2z[a(z 2 + a 2 + p 2 2 ) - 2ap 2 2 ]p 2 dp 2 dA ± 


2 2 
z c + a + 


- 


4a‘ 


2 . 2 


n?/ 2 


(A5) 


The configuration factor from the differential area dA-j^ to the 
finite area A 2 can be obtained from equation (A5) by integrating over 
P 2 , from 0 to rp which after rearranging gives 


dA->F . . „ 

^ CLA.^j 



vr 




(AS) 


The following equation, obtained by integrating equation (A 6) over 

9 1 from 0 to 2n and over z from 0 to Lq, gives the configuration 
factor from the area A^_ to the area Ap: 


ro vo vo f 
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or in nondimensional form 


Sr o' c i = 4N 


\/n 4 + 2N 2 (l + M 2 ) + (l - M 2 ) - (l - M 2 ) - N 2 


(A8) 


where 


M = r^/a 
N = l£/a 

Results obtained by evaluating equation (A8) in the range 0.1 = M ^ 1.0 
and 0.2 ^ N % 200 are given in table II and figure 4. 

A special case presents itself when becomes the full area of 

the base of the cylinder (M = l). Equation (A 8) then reduces to 




(A9) 


The reciprocity theorem can be employed to determine the configuration 
factor from the base of the cylinder to the walls F q. 

C 0' K 0 

Cones . -• When the surface of revolution is a cone (fig. 5), the 
parameters of equation (Al) can be expressed as 


cos = c - ° - g ' ^ |a - P 2 cos^G-l - e 2 ^j 


(A10) 


cos ^ 


S 


(All) 
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where 


S 2 = z 2 + p x 2 + p 2 2 - 2p x p 2 cos(9 1 - 0 2 ) (A12) 

iAg = p 2 dp 2 d0 2 

\|f = Half apex angle 
a = + z tan ^ 

Integration over 9p after substituting liquations (A10), (All), 
and (A12) into equation (Al) gives the configuration factor from the 
area dA-j_ to the differential ring in the base as 

2z cos dA-^a^z 2 + p^ 2 + p.- 2 ) - ^Pg 2 ] P 2 d Pg 


dA- F ja _ , 

1 dA.^ 2np 2 dp 2 


^z 2 + p x 2 + p 2 2 j - 4 Pi 2 p 2 2 


3/ 2 


(A13) 


The configuration factor from the differential area dA-^ to the area C-^ 
is obtained by integrating equation (A13) over p 2 from 0 to r-^. It is 


dA, F, 


~ ap i^ + p 1 2 ( 2z2 " r i 2 ) + ap -|( r -| 2 " 

2z cos i dA, J i 1 j 1 ' 1\ 1 

z 2 ) + z 2 (z 2 + r 2 ') 


r / — _ . 2 » 

-1/2 

L 4p i 

[(z > Pl 2 ) + - p^ 

1 * 


4 _ — 3 ^ o_2„ 2 _ 0 „2„ ^ _2 


P-^ - ap-j_ + 2z - az £i p 1 + z' 

4p x z 2 ^z 2 - p^j 


(Alt) 


Since the area A-^ is a surface of revolution, where 

dAf = p^ sec i|r d6 x dz 

the configuration factor from an area Rq to an area C-^ after inte- 
grating over 0j_ and collecting terms becomes 

,Lq ( a l z ^ + + c i z + d^Jdz 

T7 2 


^1, ■” 

R 0 >C 1 


^ ( agz ^ + c^.z 2 + d2Z + e 2^ 


13 2 

pliQ a^z + bjZ + c^z + d^ 

I 5 

0 a^z + b 4 z + c L 


dz 


(A15) 


ro vo MD tr* 
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where 

il 

= sec \|r 

b-^ = -3a tan \|r sec 2 4f 
°1 = (3 a ^ " r^2)tan 2 f + a 2 - r^ 2 
d 1 = -a tan t(a 2 - 

a 2 = sec \|r 

bp = -4a tan ^ sec 2 ^ 

C 2 = 2a^(2 tan 2 iV + sec^t) + 2r-^ 2 (l - tan^t) 
dp = -4a tan f(a 2 - r i^) 

e 2 = ( a2 “ r l 2 ) 2 

4 , 

= sec i|r 

p 

bj = -3a tan ^ sec ^ 
c 5 = a 2 (3 tan 2 i|c + l) 

= -a^tan \|r 

p 

a^ = see \|r 
b^ = -2a tan \|r 
2 

= a 

For brevity, equation (A15) will be written as 

F 1 
R 0’ C 1 

The first integral of equation (A15) is not an elementary integral 
and it is necessary to introduce elliptic functions in order to evaluate 
it. Reference 4 presents a method allowing integrals containing the 
square root of a quartic in the denominator of the integrand to be put 
into Legendre ' s standard form of an elliptic integral. 
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By making use of reference 5 to evaluate the elliptic integrals 
the first integral of equation (A15) becomes: 


(q - p) 

sec^k 


_ 2 («9 - ^10) . 2 (™9 - ™io)l . ^ 

I) 2 k|(x + 1) K^x 



+ 



+ E^k) + dnepgcscpg 




m=l 
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where p and q are roots of the equation 

(a + a - 3 - p)T] 2 + 2(0p - aa)ri + aa(|3 + 0) - 00(a + a) = 0 
and where a and p are the nonconjugate roots of the quartic 


z 1 * + 


S z 5 “2 -2 

io 


z<; + ^ z + ar = 0 


®2 


and 

A . sin' 1 '/ ^ 

v% + % 

CS CD = cot ® 
n n 

dn <P n Jacobi elliptic function, Vl - Ain 2 ^ 

E (jp n ,, k) incomplete elliptic integral of the second kind 

k) incomplete elliptic integral of the first kind 

G(x) = (Kcj + K^x)(Ky + Kqx) 

K 

k modulus, - 

k' complementary modulus, \Jl - k 2 

K, E complete elliptic integrals of the first and second kinds, 

respectively 

K' complete elliptic integral of the first kind with a modulus 

of k' 

Kq = a^q^ + b-^q 2 + c^q + d-^ 

K 2 = 3aqPq 2 + b 1 q 2 + 2b 1 pq + 2c 1 q + c-^p + 3d 1 
Kj = 3& 1 p 2 q + 2b 1 pq + b^ 2 + c^q + 2c x p + 3d 1 


6^ ~ Vg 

* 6*7 
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k 4 = a l p5 + b l p2 + c i p + d l 

2 

= p 2 - 2pRe(a) + | a! 

= q 2 - 2qRe(a) + J aj 

2 

Krf = p 2 - 2pRe ( p ) + | 3 ] 

Kq = q 2 - 2qRe(p) + | (3 
K 9 = -2(K^ + Kg)(Ky + Kg) 

K 10 = 5(K 5 K 8 + + 2^9) 

K 11 = " (K5K q + K^K 7 + 6K5K6) 

k 12 = 2K 6 K 8 



U = - 3% 

n 2k 


V = - 2 K 2 + 3K X 


Z(A,k) Jacobi Zeta function, E(A,k) - | F(A,k) 


^ ONONOJ 
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^2 = tan -1 


t 1 

L 0 ~ P 

£ - Lq 



OD 


irF(A,k) 
_ 2K 


The second integral of equation (A15) is an elementary integral 
which upon integration becomes 

I 2 = L 1 sin 2 t + I*) (A18) 


Equation (Al6) represents the fraction of the total radiant energy 
vhich leaves surface A x , the frustum of a right cone, and is inter- 
cepted by surface a plane circular area in the base of the cone. 

If the intercepting area A 2 is the entire base of the cone, the results 
are greatly simplified. The procedure remains unchanged up to the inte- 
gration of equation (A15) which now contains two elementary integrals 
whose integration yields 


F 

R 


1 Q 

O' 0 



csc 2 ^ 


+ 4-HL 


H 


+ CSC 


i|f(2H sin 2 y - Lq ; 


(A19) 


The configuration factor for the complete cone exchanging radiant 
energy with the complete base becomes 

F r H P = sin t (A20) 

0 

The configuration factor from the base of the cone to the walls 
can be found from the reciprocity theorem with the aid of configuration- 
factor algebra. (See appendix B. ) 

Hemisphere . - When the surface of revolution is a hemisphere (fig. 7), 
the parameters of equation (Al) can be expressed as 

g[a g - p lPg 003(6, - 8 2 )j 


aS 


(A21) 
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cos t 2 = | 

(A22) 

,2 - 2p^p 2 cos^Oj - 0 2 ^ 

(A23) 


and 


dA-^ = a d9^ dz 


dAg = p 2 dp 2 d0 2 

Integration over 02 after substituting equations (A2l), (A22), 

and (A23) into equation (Al) gives the configuration factor from the 
area dA-^ to the differential ring in the base as 


dA l F dA 1 , 2rtp 2 dp 2 


2z dA ] _ 

,( p l 2 + z2 ) 

' + ( 

' P 
z~ - 

■Pi 2 ) 

ip 2 2 ; 

OJ 

a 

r) 

OJ 

CL 

a p 2 ^ + 2p 2 2 (z 2 - 

pi 2 ) 

1 + 1 

(pi 2 ♦ * 2 ; 

f 

5/2 


(A24) 


where a is the radius of the hemisphere. Integration of equation (A24) 
over p 2 from 0 to r gives, after some rearranging, the configura- 
tion factor from the differential area dA-^ to the area C-^ as 


zr n 


dA-i F, . _ 

dA-j^, C-l 


-dA x 


i r^ + 2r^(2z 2 - e 2 ) + a 1 * 


1/2 


(A25) 


Integrating over 9-^ from 0 to 2n gives 


dA n F, 


cz 


1 2rtadz, C-^ 


(a 2 - r^ 2 ) + 


1/2 


Integrating over z from 0 to 1^ gives the configuration factor 
Rq to area as 


from area 


on o\ oj 
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L 


R 0> C l 


1 

\fkr 2 (iA 

f + ( a 2 - r 2 f . ( a 2 . r 2) 

4al^ 

V 1 \ 0 

} \ 1 / V 1 / 


(A26) 


where r-^ is the radius of the circular area in the equatorial plane, 
and Lq is the vertical height of the area R^. 

Equation (A26) may be put in nondimensional form by letting 

M = a 

» - L j/ a 


so that the configuration factor becomes 


F 

R 


Q> 1 



+ Al^N 2 - 


(l - y£)\ 


(A27) 


An interesting and useful result can be obtained from equation (A25). 
If A2 is taken as the total area of the base (i.e., r^ = a), equa- 
tion (A25) reduces to 


dA n 


dA l F cLA C 


(A28) 


By integrating over A^ the following result is obtained: 


- A P 

A l> b 0 


1 

2 


(A29) 


Now since A-j_ has not been specified it follows that the configuration 

factor from any area on the surface of a hemisphere to the equatorial 
plane is one -half. 


x 1 

Configuration Factors for the Geometry R^, Rq for the Hemisphere 

When the surfaces exchanging radiant energy are rings on the surface 
of a hemisphere (fig. 9), the parameters of equation (Al) become 



22 


cos 4r-j_ = cos ^2 


= ||l - sin 6 X sin © 2 cos^cp 1 - cpg) - cos 0 1 cos e 2 j 


(A30) 


and 


S 2 = 2 a 2 |l - sin 0i sin 0 g cos - <3g) - cos Q 1 cos 0gj (A3l) 

where 9 and q are as shown in figure 9 * 

Substitution of equations (A30) and (A 31 ) into equation (Al) yields 


dA,F 


dA^ cLAg 


1 dA x , dAg hita 2 


(A52) 


Since A 2 is a surface of revolution 

dAg = a 2 sin 0 2 d © 2 dcpg 

integration of equation (A 3 l) is simplified and gives the configuration 
factor from the differential area cLA-^ to the ring R Q as 


dA n F, 


- 


l'dA^Bj " 2 a 


(A33) 


Integration over where 

0 

dA-^ = a sin 6 -^ dG-^ 

gives the configuration factor from a ring to a ring as 

* 2**3 R 1 = rtI t > L 2 
R 2 ,R 0 


(A34) 


The areas can be expressed as 


r| = 2iral| 


Rj = 2rtaLQ 


and 
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so that 


L 1 

? = i -2 

p3 pi 2 ^ 

R 2' R 0 




P ,ii 

pl p3 2a 
2 




(A35) 


1 3 

where 1^ and are the vertical heights of the hemispherical 

segments. 
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APPENDIX B 

CONFIGURATION- FACTOR AIjGEBRA 


It is possible to determine configuration factors for many cases 
from existing configuration factors by utilizing the technique of 
configuration-factor algebra (see the section entitled "Geometric Flux 
Algebra" in ref. l) and the previously mentioned reciprocity theorem. 
The basic requirement involved in this technique is that the unknown 
configuration factors be of such a nature that they can be expressed as 
sums and differences of known configuration factors. The procedure is 
best explained by an example. 


Example 

Suppose that it is necessary to determine the configuration factor 

p 

from a ring on the surface of a cylinder R^ to an annular region 2 

of the base as shown in figure 10. It will be shown that it is possible 
to determine this configuration factor from uhe configuration factors 
for other geometries of figure 10 which are readily obtained from 
equation (A7). 

Equation (A7) gives the configuration factor for a particular 
geometry. That is, the surface of the cylinder must extend from the 
intersection of the base plane to any height above this intersection. 

The intercepting area is also restricted in uhat it must be a circular 
area the center of which is on the axis of revolution. It is therefore 
necessary to express the desired configuration factor in terms of those 
which have been obtained. This can be done as follows. From figure 10, 
it can be seen that 


and 


R i F 2 
R i> °i 


= RqF 2 

R 0' C 1 


" V 1 
° 


(Bl) 


R 1 F R 2 C = R 0 F r 2 c -¥ r i c 
1> °2 K 0> °2 O' 2 


(B2) 


1,2 " c l - 


Since C 


C 2 > then 
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R 1 F 2 

E 1 > C 1,2 


2 

= R 1 F p 
1 n 

K 1 ,0 1 


R ! F B 2 C 


(B3) 


The left-hand side of equation (B3) can be expressed as the difference 
of the right-hand sides of equations (Bl) and (B2) as follows: 


n 

K l' 2 




R 0’°2 


- R, 


0 


C R o’ c i ' Fr o' c 2 ) 


(Bl*) 


If equation (A7) is used to determine the configuration factors within 
the brackets, F,,p n can be expressed in terms of the cylinder dimen- 

R r i, 2 

sions as 


F 

R 


1 » 1 , 2 




) 


2 



Although figure 10 and the preceding example involve a cylinder, 
the procedure applies equally as well when the surface of revolution is 
a cone or a hemisphere. 

The above example is a simple application of configuration-factor 
algebra presented in order to introduce the basic ideas which will now 
be used to obtain the configuration factors for the geometries of 
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2 1 2 2 2 
R-^, C-^, R 0 ,Ci g, R-^, and R^, ^ or c yl- n< i ers > cones, and. hemi- 

5 1 

spheres as well as R;>, Rq ^ or c y llnclers 311(1 cones. (See fig. 2 for 
examples of these geometries.) 

2 

Configuration Factors for the Geometry R^, 

2 

The configuration factor for the geometry of R-pC-^ can be obtained 

from the equation derived for the geometry of R^, C^ through 
configuration-factor algebra and is given in general terms by 


R? 0 R 0’ c i 


" R n F 

0 Vl-- n 

V 


where Rq, Rq, and R-^ are ring areas. 

Cylinder . - For a cylindrical surface of revolution, equation (A7) 
can be used to determine the terms on the right-hand side of equation (B6) 
so that the configuration factor in terms of the dimensions of the cylin- 
der becomes 



Cone . - For a conical surface of revolution, equation (Al6) can be 
used to determine the terms on the right-hand side of equation (B6) so 
that the configuration factor in terms of the c imensions of the cone 
becomes 
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where the integrals 1^ and I 2 are given by equations (A17) and 
(Al8), respectively. 

The first two integrals of equation (b 8) are to be evaluated for 
a frustum of height Lq and the second two are to be evaluated for a 

frustum of height while all four integrals are to be evaluated for 

an intercepting area of radius r^. 


Hemisphere . - For a hemispherical surface of revolution, equa- 
tion (A26) can be used to determine the terms on the right-hand side of 
equation (B6) so that the configuration factor in terms of the dimensions 
of the hemisphere becomes 


F 

R 


1 ' 1 



Configuration Factors for the Geometry Rq , 2 

The configuration factor for the geometry of Rq, C^ ^ can be 

obtained from the equation derived for the geometry of Rq, ^ through 
configuration-factor algebra and is given in general terms by 


E 0’V 


= F - F 

R 0' C 1 R 0’°2 


(BIO) 


Cylinder .- For a cylindrical surface of revolution, equation (A7) 
can be used to determine the terms on the right-hand side of equa- 
tion (BIO) so that the configuration factor in terms of the dimensions 
of the cylinder becomes 


F 1 

R 0' C l, 2 





(Bll) 



28 


Cone . - For a conical surface of revoluti on, equation (Al6) can be 
used to determine the terms on the right-hanc. side of equation (BIO) 
so that the conf iguration factor in terms of the dimensions of the cone 
becomes 



Note that two integrals are evaluated for p = r , two for p = r p , 

1 " L 
and all four for a height 1,5. 


Hemisphere . - For a hemispherical surface of revolution, equa- 
tion (A26^ can be used to determine the terms on the right-hand side 
of equation (BIO) so that the configuration factor in terms of the 
dimensions of the hemisphere becomes 




p 

Configuration Factors for the Geometry of R^, C-^ ^ 

o 

Ihe configuration factor for the geometry of 2 can ^ 

obtained from the equation derived for the geometry of Rq, C-^ through 

configuration-factor algebra. This has been performed in the section 
entitled "Example" and is given by equation (B 4 ). 

Cylinder . - Equation (B 5 ) gives the configuration factor for the 
2 

geometry of R-^C-^ 2 for a cylindrical surface of revolution. 

Cone . - For a conical surface of revolution, equation (Al6) can 
be used to determine the terms on the right-hand side of equation (b 4) 
so that the configuration factor in terms of the dimension of the cone 
becomes 


ro mo vo t -1 
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Again it is to be noted that the sum of the integrals must be evaluated 
for the correct combinations of p and L. 

Hemisphere . - For a hemispherical surface of revolution, equa- 
tion (A26) can be used to determine the terms on the right-hand side of 
equation (B4) so that the configuration factor in terms of the dimen- 
sions of the hemisphere becomes 



Configuration Factors for the Geometry of 


4’ r I 


3 1 

The configuration factor for the geometry of R^Rq can 

obtained from the equation derived for the geometry of R^, C^ through 
configuration-factor algebra and is given in general terms by 


V R 1 ' i ( R 1 P r 3 pi ' R 1 P r2 c 1 

R 2' R 0 V R 1' C 0 R 1' C 0 


■ 


- R 0 F 2 

° % C 0 


(Bl6) 


Cylinder . - For a cylindrical surface of revolution, equation (A9) 
(a special case of (A7) when the intercepting area is the full base of 
the cylinder) can be used to determine the terms on the right-hand side 
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of equation (Bl6) so that the configuration factor in terms of the dimen- 
sions of the cylinder becomes 


r 3 r 1 
R 2 ,R 0 


4al| 


2I^l| + L i^( L i) + ^a 2 - Li^( L i) + 4a<: 


- l; 




4a" 


(BIT) 


Cone . - For a conical surface of revolution, equation (A19) (a 
special case of equation (Al6) when the intercepting area is the full 
area of the base of the cone) can be used to determine the terms on the 
right-hand side of equation (Bl6) so that the :onf iguration factor in 
terms of the dimensions of the cone becomes 


R 2' R 0 2l| (l R + L R ) _ 




2 , H H 

esc t + 4L^L^. - L 


i\/( I if csc2 * + 


- Lq \j(4>) cs °4 + ^ hl 3 + 4f4) cs ;2 + + 


+ 2 


COS ^ cot \|M 


(B18) 


Hemisphere . - For a hemispherical surface of revolution, equa- 
tion (^A2^5 cannot be applied to equation (Bl6) to determine the config- 

3 1 

uration factor for the geometry of This is because the first 

two terms on the right of equation (Bl6) cannot be obtained from equa- 
tion (A26), which is for hemispheres, since the surface of revolution 
above plane 1 is a spherical segment. Therefore, the equation for the 
configuration factor between rings on the surface of a hemisphere is 
derived in appendix A (eq. (A35)) by integrating equation (4) over 

Rq and R^. 


(-3 o\ CT\ cu 
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O p 

Configuration Factors for the Geometry R^, R 1 

Radiation in a closed system composed of several isothermal surfaces 
is further complicated by surfaces of positive curvature. A surface of 
positive curvature Aj intercepts a portion of its own radiated energy 

and a configuration factor of the form must be determined. In 

order to determine F^j, use is made of the fact that if a surface Aj 

is surrounded by (n - l) other surfaces, forming a closed system, the 
sum of the configuration factors between Aj and all surfaces is one, 
or 


I 

1r=1 




= 1 


(B19) 


If the surface Aj can see itself, equation (B19) can be written as 


j-1 


n 


f j, j ~ 1 






(B20) 


k=l k=j+l 

For a surface of revolution, equation (B20) can be replaced by 


Fp p - 1 - F p i - F? p 
R^.RV R, ,ci R ^ C o 


(B21) 


'l' x l 


1’ ^0 


P 

where R^ is the ring formed by the intersections of planes 1 and 2, 

1 2 

which are perpendicular to the axis of revolution. Areas C q and Cq 

are circular areas in planes 1 and 2 bounded by the surface of revolu- 
tion. The last two terms in equation (B2l) represent the fractions of 

2 2 

radiation from R-^ which fall on surfaces below and above R^, 

respectively. 

Cylinder . - If the surface of revolution is a cylinder 


Co 




and equation (B2l) reduces to 


? 2 2 

R i' R i 


2F 


f *' 1 




(B22) 
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The configuration factor F 2 1 from the walls of a cylinder to 

R l' C 0 

the base can be obtained from equation (A 9 ), and equation (B22) may then 
be written in terms of the dimensions of the cylinder as 


R 2 R 2 
R 1' R 1 


2 a 


1 + 


1 2 a y 


(B23) 


It can be seen from equation (B23) that F 0 n for a cylinder is 

R 1' R 1 

dependent upon the height of the ring but not upon the position of the 
ring above the base of the cylinder. 


Cone . - For a conical surface of revolutior 

F 


R 2 r 2 
R 1' C 0 




(B24) 


where is all the surface area above R^. 

By using equation (B24), equation (B2l) becomes 


F r 2 r 2 " 1 ' F r 2 C 1 " F £ pH 


(B25) 


The portion of the energy leaving R-^ which is intercepted by the 
circular area Cq can be determined from equation (A19) as 


r p 1 “ 
R l’°0 2 


$♦*§) 


csc 2 \|r + + esc i^ 2 I^ sin 2 \|/ - L 2 


(B26) 


Configuration- factor algebra for the exchange of radiant energy 
between and R§ gives 


R 1 F 2 it = R 2 F h 2 = R 2 F H 2 ' K F JE 1 “ R 1 F 2 1^ 

1 Rf,R| 2 R R Ri l|,C§ l Hj,cJ R l> C 0, 


(B27) 
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By using equations (A 19 ) and (A20), equation (B 27 ) yields 


2(1? ♦ if 


^fcsc 2 * + 4L R I^ - CSC 4^2L^sin^V + L^j 


(B28) 


Substitution of equations (B26) and (B28) into equation (B25) gives 
the configuration factor for a ring of a cone radiating to itself as 


F =1 

? ? x 

R I' R 1 


t H t H 


ur csc^ cos ^ cot ^ 


(B29) 


It can be seen from equation (B29) that the portion of its own 
radiated energy that a ring on the surface of a cone receives is dependent 
upon the relative position of the ring as well as on the height of the 
ring. 


Hemisphere . - For a hemispherical surface of revolution, since the 
total height equals a, equation (B 2 l) can be written 


F P p = 1 - F 0 - F Q d 

r !’ r ! R r c J R ?’ R I 


(B30) 


where 


P 

RqF 


R 2 C 1 
K l> c 0 


R 0’ C 0 


"o* 1 + R 1 R 2 „1 

R 0> c o R i' R o 


(B31) 


By using equations (A29) and (A 35 ), equation (B30) becomes 


R 2 ,R 2 " 2 a 


(B32) 


Here again, as for a cylinder, F 2 2 ^P^^ent upon the 


R i' R l 


height of the ring and not upon its position above the base plane. 



TABLE II.- CONFIGURATION FACTORS FOR GEOMETRY Ri,C n FOR CYLINDERS 
















TABLE III.- CONFIGURATION FACTORS FOR GEOMETRY R^C-. FOR CONES 
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